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Abstract 

We study the emission and absorption spectra due to various photon 
and pair processes in a non-equilibrium pair plasma containing a signifi- 
cant density of photons. We present here some preliminary results from 
Monte-Carlo simulations. These investigations are likely to be useful in un- 
derstanding the radiation and relaxation mechanisms in non-thermal 7-ray 
sources in astrophysics. 
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1 Introduction 



Relativistic plasmas composed of electron-positron pairs, ions, and photons are 
believed to arise from various astrophysical processes involving neutron stars and 
black holes. There is some indirect evidence for the existence of such plasmas 
deduced from observations of active galactic nuclei and various 7-ray sources. 
Such plasmas have been studied extensively by many authors in the past. See for 
example [1 — 3] and references therein. Much of the past research was devoted 
to the study of equilibrium plasmas with internally generated photons. It will be 
interesting to study the evolution of a relativistic non-thermal pair plasma with 
an initial distribution of high-energy (0.1 — 100 MeV) photons, whose density is 
comparable to the pair density, which are not produced by any emission processes 
within the plasma. This will be of some use in imderstanding the physics of non- 
thermal 7-ray sources in astrophysics. It may also have some relevance to the 
study of 7-ray-burst sources in their final stages when they are only moderately 
optically thick. The present work is an attempt to extend and generalize the 
existing research on equilibrium pair plasmas mentioned above. 

2 The Model 

We consider a non-equilibrium plasma composed of photons, protons, and rel- 
ativistic electrons and positrons. The plasma is assumed to be stationary (i.e., 
confined) and free from any magnetic fields. The number densities of photons, 
electrons, and positrons are n^, n_, and n+, respectively. They are assumed 
to be comparable to each other. In such a case many radiation and relaxation 
mechanisms have to be taken into account to determine the evolution of the 
plasma. Important photon processes are Compton (and inverse Compton) scat- 
tering, double Compton scattering, and pair creation. The particle processes 
are bremsstrahlung from e^-e^, e^-e^, and electron-proton collisions, M0ller and 
Bhabha scattering, and pair annihilation. 

2.1 The Kinetic Equation 

To study the evolution of such a plasma we proceed from the relativistic Boltz- 
mann equation, 

p''^,f{x,p)=Y.[71,-X^f], (1) 
i 

where and are the position and momentum 4-vectors, and / is the dis- 
tribution function. The sum is over all processes that contribute to emission 
and absorption, and rji and Xi ^^e the corresponding binary coeffiecients. For 
isotropic emission and absorption, we define the spectral functions for photons 
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and particles, 

4 7J- ^2 4 TT 7 - 1) 

F,{e) = ^f,{e) and ^±(7) = (2) 

respectively. Here e is the photon energy in units of mc^ and 7 is the particle 
Lorentz factor. They are normalized so that 

J deF^{e) = l and J d^F±{^) = l. (3) 

2.2 A Family of Non-Equilibrium Functions 

We use the following three-parameter family of non-equilibrium spectral func- 
tions: 

= i '""P^^/"^ (4) 

C2e^^, otherwise. 

Here a £max represents the spectral break, where £max is the energy at which the 
thermal part has its maximum, and G = kBT/mc^ is the "temperature" of the 
low-energy part of the spectrum, where m is the electron rest mass and A;^ is 
the Boltzmann constant. Finally, S is the power-law index for the high-energy 
part. The constants ci and C2 are determined by the continuity of the function at 
a £max and the normalization given by equation (3) . A similar family of functions 
is defined for particles as well. 

The motivation for this choice is the following: It has been realized for some 
time'' that M0ller and Bhabha scattering processes are not efficient enough to 
maintain particle equilibrium in a relativistic plasma for > 3.5. This result 
applies to pair-dominated plasmas. However, if we also include energetic photons 
of comparable density which are non-thermal, then it is likely that a non-thermal 
tail will develop beyond a few MeV for pairs as well as photons. We parametrize 
this by a power-law index S. Besides, similar power-law tails are found in many 
7-ray-burst spectra'^. The low-energy part is characterized by a fictitious tem- 
perature 0. This is because at low energies many of the relaxation processes are 
efficient for thermalization. This is only a hypothetical family of functions; we 
do not imply that they are based on any confirmed observations. 

2.3 The Absorption and Emission Coefficients 

In binary collisions between two species (1 and 2) of particles or photons with 
spectral functions Fi and F2, the photon emission coefficient is given by, 

^ ^ = 1« 2M ^ / 11 (^^^i ^^^^^ ' ^ 

167r^(l + 0i2j JR jt=i ^1^2 dF' 



3 



Here dP = s'^dedfl for photons and 'y'^pdjdfl for particles, dai2/dP is the 
differential cross section for the process, 5i2 = 1 for identical particles and is 
zero otherwise, c j3rei is the relative velocity of the particles, e-i = e in the case of 
photons and = 7 in the case of particles, and pi are the 4-momenta of particles 
or photons in units of mc. The integration is over all initial states, subject to 
the conservation of 4-momentum so that one of the final states is a photon of 
energy e as required. This usually results in a multi-dimensional integral over 
a complicated region R. The photon absorption coefficient is obtained from a 
similar formula. In that case the integration is performed over all final states 
which absorb an initial state of energy e. Physically, 47r£rj/mc gives the rate 
at which the photons are emitted per unit volume, per unit time, and per unit 
energy. An analogous interpretation is given to the photon absorption coefficient 
and to the particles cases. 



3 Results and Discussion 

So far we have performed Monte-Carlo evaluations of various emissivity integrals 
for pair annihilation, Compton scattering, and e^-e^ and e^-e^ bremsstrahlung. 
To test our calculations we have used relativistic Maxwell-Boltzmann distribu- 
tions at various temperatures and compared our results on equilibrium plasmas 
with those of several previous authors^^^; our results are found to be in good 
agreement. Using the same programs, we have calculated the spectra for the 
family of non-equilibrium functions specified by equation (4). We define a di- 
mensionless emissivity function. 

Til n2 Tg 

where rg is the classical electron radius. The function \I/ due to pair annihilation 
and Compton scattering is shown in Figure 1. We are presently computing the 
emission and absorption functions for the remaining photon and particle inter- 
actions using the same Monte-Carlo scheme. This will enable us to understand 
various relaxation and radiation mechanisms in such plasmas. We will report 
these investigations in greater detail elsewhere^". 
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Figure 1: Function for (a) pair annihilation and (b) Compton scattering. In 
all cases we use a = 3 and ^ = 2, 6 = 0.1, 0.33, 1, 3.3 and 10. 
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